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Abstract. We construct an expanding gradient Ricci soliton in dimension three over the topo- 
logical manifold RxT 2 that aproaches asymptotically a constant curvature cusp at one end, and 
a flat manifold on the other end. We prove that this is the only gradient soliton with this topol- 
ogy, provided the curvature is negatively pinched, —1/4 < sec < 0, at the time-zero manifold 
(normalizing the soliton to be born at time —1). 



1. Introduction 

A gradient Ricci soliton is a smooth riemannian metric g on a manifold M together with a 
potential function / : M — > M. such that 

(1) Ric + Hess/+ |g = 

for some e G R. Solitons provide special examples of self-similar solutions of Ricci flow, Jj<7 = 
— 2Ric, evolving by homotheties and diffeomorphisms generated by the flow (/>(■, t) of the vector 
field grad/, this is g(t) = (et + l)(j>*go. The constant e can be normalized to be — 1, 0, 1 according 
the soliton being shrinking, steady or expandig respectively (see [4] for a general reference) . Solitons 
play an important role in the classification of singular models for Ricci flow despite of (or actually 
due to) existing only a limited number of examples. In dimension 3, the only closed gradient 
solitons are those of constant curvature. Furthermore, by the results of Hamiton-Ivey [S], [TT] and 
Perelman [13] , the only three-dimensional open gradient shrinking solitons with bounded curvature 
are S 3 , R 3 and § 2 xl with their standard metrics, and their quotients. Notice that in all these 
examples the gradient vector field is null. Examples with nontrivial potential function in dimension 
three include the Gaussian flat soliton, the steady Bryant soliton, the product of the 2-dimensional 
steady Cigar soliton with R due to Hamilton, and a continuous family of rotationally symmetric 
expanding gradient solitons due to Bryant. In summary, shrinking and steady solitons are very 
few, and these are useful in the analysis of high curvature regions of the Ricci flow. Expanding 
solitons are less understood and there is much more variety of them. 

A couple of motivating examples are the following. The hyperbolic metric g = dr 2 + e~ 2r (dx 2 + 
dy 2 ) on R 3 together with a trivial potential / = est fits into the soliton equation with e = 1, so 
any quotient (hyperbolic manifold) yields an expanding soliton. An open quotient of the hyperbolic 
space, for instance the cusp RxT 2 obtained as a quotient by parabolic isometries (represented by 
euclidean translations on the xy-plane) yields also a soliton. 

We will restrict ourselves to bounded curvature metrics that yield uniformly bounded curvature 
flows. From the PDE point of view, this condition ensures existence and uniqueness of solutions 
of the flow, both in the compact case ([8], [6]) as well as on the noncompact one ([15], [3j). If 



Universitat Autonoma de Barcelona, Departament de Matematiques - 80193 Bellaterra, Barcelona (Spain)- 
E-mail: dramosSmat . uab . cat . 



2 



DANIEL RAMOS 



the uniformly boundedness condition of the curvature is dropped, we can loose the uniqueness; 
for instance approximating a cusp by high-curvature capped ends (cf. with 16 and [14 in the 2- 
dimensional case). Nevertheless, even with this assumption, it is not clear what an open expanding 
hyperbolic manifold is at the birth time, namely t — > —1 when the evolution is g(t) = (t + l)gh yp - 
A sequence of shrinked negatively curved manifolds does not need to have a limit in the Gromov- 
Hausdorff sense, since the curvature is not uniformly bounded below. A hyperbolic cusp on 1 x T 2 
as an expanding soliton tends to a line while its curvature tends to — oo, as t — > —1 the birth time. 

Another interesting phenomenon occurs in K™ endowed with the euclidean metric: it fits into 
the soliton equation together with a potential function f(p) = —£4- for any e € R. The nonzero 
cases are the so called Gaussian solitons, and even when the metric is constant in all cases (hence 
there is a unique solution with a given initial condition), there is more than one soliton structure 
on it. 

The aim of this paper is constructing a particular example of expanding gradient Ricci soliton on 
RxT 2 , different from the constant curvature examples. Furthermore, we prove that it is the only 
possible nonhomogeneous soliton on this manifold provided there is a lower sectional curvature 
bound equal to — \. 

In Section [2] we consider the generic metric g — dr 2 + e 2h (dx 2 + dy 2 ) over M = R x T 2 , where 
h = h(r) is a function determining the size of the foliating flat tori, and a potential function 
/ = f( r ) constant over these tori. We find a suitable choice of h and / that makes the triple 
(M, g, f) a soliton solution for the Ricci flow with bounded curvature, by means of the phase 
portrait analisys of the soliton ODEs. 

Theorem 1.1. There exists an expanding gradient Ricci soliton (M,g,f) over the topological 
manifold M = R X T 2 satisfying the following properties: 

(1) The metric has pinched sectional curvature — \< sec < 0. 

(2) The soliton approaches the hyperbolic cusp expanding soliton on one end. 

(3) The soliton approaches locally the flat Gaussian expanding soliton on a cone on the other 
end. 

More preciselly, M admits a metric 

g = dr 2 + e 2h ^(dx 2 + dy 2 ) 

where (r, x, y) £ R x S 1 x S ; and a potential function f = f(r), satisfying the soliton equation and 
with the stated bounds on the curvature, such that 

r 



and 



and / — > est as 



h ~ In r and f ~ as 

J 4 



For the asymptotical notation we write 

4>{r) ~ tp(r) as r — > 00 

if 

r->oo 1p(r) 

Let us remark that when r — > +00 the theorem states that the metric approaches g = dr 2 
r 2 (dx 2 + dy 2 ). This is a nonflat cone over the torus, namely its curvatures are sec rx = sec ry = 
and sec xy = —^s, but it indeed approaches a flat metric when r — > +00. 

This example is interesting in regard of the following known fact. 
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Proposition 1. Let (M n ,g, f) with n > 3 be a complete noncompact gradient expanding soliton 
with Ric > ^ + 5 for some 5 > 0. Then f is a strictly concave exhaustion function, that achieves 
one maximum, and the underlying manifold M n is diffeomorphic to R™. 

Cf. with [SJ, Lemma 5.5 and Remark 5.6. Our example proves that the lemma fails if one only 
assumes Ric > In this critical case the soliton also has a strictly concave potential, but in 
this case / has no maximum and it is not exhausting, and actually this solution admits a different 
topology for the manifold, namely RxT 2 . 

In Section [3] we consider the general case of a metric over IxT 2 with sec > — \ , and we 
prove that the only nonfiat solution is the example previously constructed. The lower bound on 
the curvature implies a concavity property for the potential function. This leads together with 
the prescribed topology to a general form of the coordinate expression of the metric, that can be 
subsequently computed as the example. 

Theorem 1.2. Let [M,g,f) be a nonfiat gradient Ricci soliton over the topological manifold 
M = M x T 2 with bounded curvature sec > — t ■ Then it is the expanding gradient soliton depicted 
on Theorem (1.1). 



In section [4] we explore a growth property of the scalar curvature on our soliton. In nonnegative 
sectional curvature, the evolution of Rvn and Harnack inequalities [TO] imply that > — y 
pointwise, and in particular ^ > for all t if the solution is also ancient. Our example exposes 
that this is not the case in negative scalar curvature, even with a soliton solution. Despite the 
self-similarity, the behaviour of the curvature growth is different at different times. The combined 
effects of the diffeomorphism translation and the homothety act in opposite manner. For short 
time after birth, the negative curvature is increasing everywhere. After some small time, there 
appear points where the curvature is decreasing, but eventually all points recover the increase of 
the scalar curvature and the limit of the curvature is zero for every fixed point as t — > +00. 

Theorem 1.3. Let (M,g(t)) be the (soliton) Ricci flow defined on M = 1 x T 2 and for t G 
(— l,+oo), such that g(0) = go where go is the metric constructed in Theorem \l.l\ Let R = R(t) 
be the scalar curvature of g(t). Then there exists S < such that 

• for all t G (5, +00) there exist points in M with ^ > and points with ^ < 

• for some —l<t<5, it is satisfied ^? > everywhere in M. 



Most tedious computations thorough the paper can be performed and checked using Maple or 
other similar software, therefore no step-by-step computations will be shown. Pictures were drawn 
using Maple and the P4 program pQ. 

Acknowledgements: The author was partially supported by Feder/Mineco through the Grant 
MTM2009-0759. He is indebt to his advisor, Joan Porti, for all his guidance. He also wishes to 
thank Joan Torregrosa for pointing out the technique of compactified phase portraits. 

2. The asymptotically cusped soliton 
Let us consider the metric 
(2) g = dr 2 + e 2h(r \dx 2 + dy 2 ) 

where h = h{r) is a one-variable real function, and a potential function / = f(r) depending also 
only on the r-coordinate. The underlying topological manifold can be taken (r,x,y) G K x S 1 x S 1 = 
RxT 2 since M 3 with this metric admits the appropriate quotient on the x,y variables. Standard 
riemannian computations yield the following equalities. 
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Lemma 1. The metric in the form associates the following geometric quantities: 

Ric = -2((h') 2 + h")dr 2 - e 2h (h" + 2{h') 2 )(dx 2 + dy 2 ), 
R = -Hi" - 6(h') 2 , 
sec xy = -(h 1 ) 2 , 
sec rx — sec ry = —{{h 1 ) 2 + h"), 

Hess / = f'dr 2 + e 2h fh'(dx 2 + dy 2 ), 

Af = 2h'f + /", 
grad/ = /'&, 
V/ - f'dr, 
|grad/| 2 = (/') 2 . 



Hence the soliton equation ([T]) for this metric turns into 

(| + /" - 2{h') 2 - 2/i") dr 2 + e 2h (| + ti f - 2(h') 2 - ti'^j {dx 2 + dy 2 ) = 0. 

This tensor equation is equivalent to the ODEs system 

f f + /" - 2(/ l ') 2 - 2h" = 
W \ | + V/' - 2{h') 2 - h" = 0. 

Let us remark that this system would be of second-order in most coordinate systems, but in ours 
we can just change variables H = h' and F = f , and rearrange to get a first-order system 

( H' = HF- 2H 2 + f 
\ F' = 2HF - 2H 2 + ^. 

We can solve qualitativelly this system using a phase portrait analysis (see Figure [T]). Every 
trajectory on the phase portrait represents a soliton, but will not have in general bounded curvature. 
Actually, bounded curvature is achieved if and only if both H and H' are bounded on the trajectory. 

The critical points (stationary solutions) of the system are found by solving {H' = 0, F' = 0}. 
If the soliton is shrinking (e = —1), there are no critical points and no trajectories with bounded 
curvature, agreeing with Perelman's classification. If the soliton is steady (e = 0), there is a whole 
straight line {H = 0} of fixed points representing all of them the flat steady soliton. In this 
case there are neither trajectories with bounded H, hence all solutions have unbounded negative 
curvature at least in one end. Let us assume henceforth that the soliton is expanding (e = 1), so 
our system is 

f H> = HF-2H 2 + \ 
1 ' \ F' = 2HF-2H 2 + \. 

There are two critical points, 

(H,F) = (±±0). 

The critical point (s,0) corresponds to a soliton with h(r) = | + c\ and f(r) — C2, the gradient 
vector field is null, and the metric is go = dr 2 + e r+Cl (dx 2 + dy 2 ), which is a complete hyperbolic 
metric, with constant sectional curvature equal to — |, and possesses a cusp at r — > — oo. As a 
Ricci flow it is g(t) = (t + l)go, it evolves only by homotheties, and it is born at t = — 1. The 
symmetric critical point (—^,0) represents the same soliton, just reparameterizing r — > —r. 
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Figure 1. Phase portrait of the system Q. 

The phase portrait of the system Q has a central symmetry, that is, the whole phase portrait is 
invariant under the change (H, F, r) — > (-H, —F, — r), so it is enough to analyze one critical point 
and half the trajectories. 

We shall see that the critical points are saddle points, and there is a separatrix trajectory 
emanaing from each one of them that represents the soliton metric we are looking for. Both 
trajectories represent actually the same soliton up to reparametcrization. 

Lemma 2. Besides the stationary solutions, and up to the central symmetry, there is only one 
trajectory S with bounded H. This trajectory is a separatrix joining a critical point and a point in 
the infinity on a vertical asymptote. 

Proof. The linearization of the system is 

/ H 1 \ _ f F- AH H \ f H \ 
\ F' ) ~ \ 2F-AH 2H ){ F J' 

The matrix of the linearized system has determinant —4H 2 < 0, so the critical points are saddle 
points. For each one, there are two eigenvectors determining four separatrix trajectories; being 
two of them attractive, two of them repulsive, according to the sign of the eigenvalue. 

We are interested in one of the two repulsive separatrix emanating from the critical point 
(H, F) = (A 0), pointing towards the region {H < ~,F < 0}. Wc shall see that has this is the only 
solution curve (together with its symmetrical) with bounded H along its trajectory, so it represents 
a metric with bounded curvature. 

In order to analyze the asymptotic behaviour of the trajectories, we perform a projective com- 
pactification of the plane, as explained for instance in [7], Ch. 5. The compactified plane maps into 
a disc where pairs of antipodal points on the boundary represent the asymptotic directions, Figure 
[2] shows the compactified phase portrait of Q . A standard technique for polynomial systems is to 
perform a change of charts on the projective plane so that critical points at infinity can be studied. 
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Figure 2. Compactified phase portrait of the system Q. 



A sketch is as follows: a polynomial system 



{ 



x = P{x,y) 

y = Q(x,y) 



can be thought as lying on the {z — 1} plane in the xyz-space. By a central projection this maps 
to a vector field and a phase portrait on the unit sphere, or in the projective plane after antipodal 
identification. In order to do this, it may be necessary to resize the vector field as 



so that the vector field keeps bounded norm on the equator. However, this change only reparame- 
terizes the trajectories. A global picture can be obtained by orthographic projection of the sphere 
on the equatorial disc, as in Figure [5] or it can be projected further to a plane {x — 1} or {y = 1} 
in order to study the critical points at the infinity. Let us remark that this technique works only 
for polynomial systems since the polynomial growth ratio suits the algebraic change of variables. 
In our system, this analysis yields that for every trajectory the ratio H/F tends to either 0, 

1 + Y or 1 - ^ as r -> ±oo; represented by the pairs of antipodal critical points (of type 
node) at infinity. The knowledge of the finite and infinite critical points, together with their type, 
determines qualitatively the phase portrait of Figure [2] by the Poincare-Bendixon theorem. Thus, 
a trajectory with bounded H on the M 2 portrait, when seen on the RP 2 portrait must have their 
ends either on the finite saddle points or on the infinity node with H/F ratio equal to (meaning 
a vertical asymptote). The only trajectory satisfying this condition is the claimed separatrix and 
its symmetrical. □ 

We shall see that this trajectory S is parameterized by r £ (— oo, +oo), and when r — > — oo the 
function h(r) behaves as | and then the solution is asymptotically a cusp. Similarly, we will see 
that h' , h" — > when r — > +oo and then the solution is asymptotically flat. 

To better understand the phase portrait it is useful to consider some isoclinic lines. This will 
give us the limit values for H, H' and the range of the parameter. 

Lemma 3. The vertical asymptote for the trajectory S occurs at H — 0. Furthermore, it is 
parameterized by r € (— oo, +oo) and H, H' — > as r — > +oo. 
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Proof. The vertical isocline {H' = 0} is the hyperbola 

f = 2H -Jh 

and the trajectories cross it with vertical tangent vector. The horizontal isocline {F' = 0} is the 
hyperbola 

F=U2H 1 



2H y 

and the trajectories cross it with horizontal tangent vector (see Figure |3| . An oblique isocline is 
the hyperbola 

since over this curve the vector field has constant direction: 

{H',F')\ {HAH _^ = i,6i/ 2 - 3 -) = l -) ( I . •'-> ) . 




Figure 3. Close-up of the separatrix trajectories (bold lines); the vertical isocline 
(dashes and dots); and the horizontal and oblique isoclines (dots and arrows). 

All three isoclines intersect at the critical points. Furthermore, the tangent directions at the 
critical point (^,0) have slope ^jj\ H _i = 4 for the vertical isocline, ^jj\ H _i = 2 for the horizontal 

one, and = 8 for the oblique one. The separatrix lines emanating from the critical point 

follow the directions given by the eigenvectors of the matrix of the linearized system (evaluated at 
the point), that is, the matrix 



- 2 2 



i 

2 

-2 1 

whose eigenvalues are ~ 1 ^" v/ ^ and ~ 1 ~ v/ ^ with eigenvectors 

1 A / 1 



3 + V5 J ' 1 3-V5 
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respectively, so the separatrix lines have tangent directions with slope 3 + y/E ~ 5.24 and 3 — y/5 ~ 
0.76 respectively. The repulsive separatrices are the ones associated with the positive eigenvalue, 
that is and the slope is 5.24. 

The repulsive separatrix emanating towards {H < |, F < 0} initially lies below both the vertical 
and horizontal isoclines, so it moves downwards and leftwards; and above the oblique one. The 
horizontal and oblique isoclines form two barriers for the separatrix, this is, the separatrix cannot 
cross any of them. This is obvious for the horizontal one, since the flow is rightwards and the 
trajectory is on the right. For the oblique isocline, we just check that any generic point on the 
isocline (H,AH — w) has tangent vector (1,4+ -pp) and a normal vector v — (—4— -g^, 1) pointing 
leftwards and upwards for < H < | . The scalar product of the normal vector v and the vector 
field (H' , F') over this isocline is 

= { 2H2 -\) (- 4 - h + 3 ) = - 2H2 + 2jp-l >0 

whenever < H < i. This means that the flow is always pointing to the left-hand side of 
the isocline branch and therefore is a barrier. This proves that the separatrix moves downwards 
between the two barriers and therefore H — > 0. 

Actually, the vertical isocline is also a barrier for the separatrix. For, if at some point it touched 
the vertical isocline, it would then move vertically downwards, keeping the trajectory on the right- 
hand side of the isocline. There would be then a tangency, but it is impossible since the tangent 
vector should be vertical. Since the vertical isocline is becoming itself vertical, this means that 
the vertical isocline acts as an atractor for the trajectories. Indeed, the trajectory lies initially in 
the region {H > 0, H' < 0} but H must remain positive since it cannot cross the vertical isocline. 
Therefore H is positive and decreasing, so H' must tend to 0. This implies that the trajectory 
tends to the vertical isocline. It is important to note that both the vertical and horizontal isoclines 
come close together when H — > 0, but the trajectories stick to the vertical one much faster than 
to the horizontal one. 

We now see that r G (— oo, +oo). This follows inmediately from the Hartman-Grobman theorem 
for the case r — > — oo, but the trajectory might, a priori, escape to infinity in finite time. This would 
require that the velocity tangent vector tends to infinity in finite time, but this is impossible, since 
H and H' are bounded, thus F' is bounded and hence the tangent vector (H' , F') is bounded. □ 

We have seen that not only H —> as r —> +oo, but also that H' — > 0. That is, hi , h" —> 
as r — > +oo, so all the sectional, scalar and Ricci curvatures tend to zero, the metric becoming 
asymptotically flat. 

At this point we have seen the existence of the soliton asserted in Theorem |I.1| We now give 
some more detailed information about the asymptotic behaviour of / and h at the ends of the 
manifold. 

Lemma 4. The asymptotic behaviour of f and h is 

t 



and f — > est as 
2 



and 



h In r and / ~ 
Proof. Recall a version of the l'Hopital rule: if 



4 



lim 4>(x) = lim ip(x) — 0, ±oo and lim 

x— foo x— yoo x^roo w (x) 
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then 

The case when r — > — oo follows from Hartman-Grobman theorem: the phase portrait in a 
small neighbourhood of a saddle critical point has a flow that is Holder conjugate to the flow of a 
standard linear saddle point 

x — x 



a T \ a 



y = -y, 

whith solution x(t) = fcie*, y(t) = kie~ l . This means that r is defined from — oo onwards, that 
if -> \ and F, H', F' -> and ->• 3 + \/% as r -> -oo. Since H = ti -> \ then h -> -oo 

and 1- By l'Hopital, ^ ~ ^ 1- Using more accurately the Hartman-Grobman theorem, 

there exists a Holder function rj : [7 C R — s- K defined on a neighbourhood of zero, and constants 
a, C > such that F(r) = r)(kie r ) and 

|F(r) - F(r )\ = \v(he r ) - r^e'' )! < C\k x e r - k x e r 
When rg — > — oo, we obtain 

\F(r)\ < Ce c 

thus, F is integrable on an interval (— oo, c] and thus / — > /o = csi as r -> — oo. The constant /o is 
actually sup / and can be chosen since it bears no geometric meaning. A more accurate description 
of / using l'Hopital tells 

h-r/2 r 
lim '— = 3 + V5. 

r->-°° f - / 

For the case when r — > +oo, we know from the trajectories that if, if' — » 0. Letting r — > +oo 
in the first equation of (jij we deduce that HF — > — |. Using this and letting r — > +oo in the 
second equation of (W| we conclude that F' — > — |, that is, F" ~ — | and by l'Hopital, F ~ — ~ 

2 — 

and / ~ — ^- as r — > +oo. 

Now, since lim,._>. +00 — = — |, we have 

1 rrr ^ ,. F 1 if 
= hm if F = hm — — : — = hm — — r 

2 v— *-+oo r— >+oo 7" ^ 7* 2 r— >+oo 7" -'- 

thus if ~ i and therefore /i ~ lnr as r — > +00. □ 

r 

It remains only to check the bounds on the sectional curvatures. 

Lemma 5. The metric with the function f obtained as solution of the system has bounded 
sectional curvature 

— - < sec < 0. 
4 

Proof. The expression for the sectional curvatures is given in Lemma [T] The case 

sec xy = -{h') 2 = -if 2 

is trivial since < if < \ and therefore — \ < sec xy < 0, tending to — \ on the cusp end, and to 
on the wide end. The other sectional curvatures are 



= -{{h'f + h") 



= -ff 2 - if' 

" 2 



= if 2 - HF - 1 



10 



DANIEL RAMOS 



We saw in Lemma[3]that {F 1 = 0} = {2HF-2H 2 + \ = 0} is a barrier for the separatrix S. Hence, 
F' < along S and therefore sec rx , sec ry > — \ . Similarly, the set {H 2 — HF — \ = 0} can also be 
checked to be a barrier for S (actually a barrier on the opposite side), and hence sec rx , sec ry < 0. 
We also saw in Lemma [ij the asymptotics of F', therefore we have — \ < F' < and sec rx — sec ry 
tend to — | on the cusp end, and to on the wide end. □ 

This finishes the description of the soliton stated on Theorem 

3. Uniqueness 

Let (M,g, /) be a gradient expanding Ricci soliton over R x T 2 such that sec > —1/4. Then 

Ric + Hess/ + ^g = 0, 

sec > 1/4, 
Ric > -1/2, 
R > -3/2. 

Recall a basic lemma about solitons, that can be proven just derivating, contracting and commuting 
covariant derivatives on the soliton equation, see [4]. 

Lemma 6. It is satisfied 

R + A/ + 3/2 = 0, 
#(gradi?, •) = 2Ric(grad/, •), 
i?+|grad/| 2 + / = a 

Since the soliton is defined in terms of the gradient of /, we can arbitrarily add a constant to / 
without effect. We use this to set C = —3/2 above so that we have 

A/ = |grad/| 2 + /. 

The bound on the curvature implies 

Hess / < 0, 

A/<0, 

(gradi?,grad/) > -|grad/| 2 . 

First equation means that / is a strictly concave function (— / is a strictly convex function), i.e. 
— /07 is a strictly convex real function for every (unit speed) geodesic 7. This is a strong condition, 
since then the superlevel sets A c = {/ > c} are totally convex sets, i.e. every geodesic segment 
joining two points on A c lies entirely on A c . Second equation is just a weaker convexity condition. 
This concavity on this topology implies that / has no maximum. 

Lemma 7. The function f is negative and has no maximum. 

Proof. Note that / is bounded above since / = A/ — | grad f\ 2 < 0. Now suppose by contradiction 
that the maximum of / is attained at some point of RxT 2 , then we can lift this point, the metric 
and the potential function to the universal cover 1 x I 2 . There is then a lattice of points in 
the cover where the lifted function / attains its maximum. But this is impossible since a strictly 
concave function cannot have more than one maximum (the function restricted to a geodesic 
segment joining two maxima would not be strictly concave). 

□ 
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Remark. As stated in Proposition [TJ if Ric > — | + S for any 5 > 0, then / has a maximum, and 
the set A = {/ > f max — fi} is compact and homeomorphic to a ball for small /i. The function / is 
then an exhaustion function, this is, the whole manifold retracts onto A via the flowline of / and 
therefore M = K 3 . Thus this stronger bound on the curvature is not compatible with M = E x T 2 . 

Now we prove that level sets of / are compact. 

Lemma 8. The function f is not bounded below and the level sets {/ = c} are compact. 

Proof. Consider M — T 2 x K as splitted into T 2 x (— oo, 0]UT 2 x [0, +oo), each component containing 
one of the two ends. Since / has no maximum, there is a sequence of points {x{\ tending to one 
end such that f(xt) —> sup /. Let us assume that this end is T 2 x (— oo, 0]. Then when approaching 
the opposite end / is unbounded. Indeed, suppose by contradiction that there is a sequence of 
points {yi} tending to the +oo end such that f{y,i) — > —K > — oo. There is a minimizing geodesic 
segment ji joining Xi with y±. This gives us a sequence of geodesic paths (whose length tends 
to infinity), each one crossing the central torus T 2 x {0}. Since both the torus and the space of 
directions of a point are compact, there is a converging subsequence of crossing points together 
with direction vectors that determine a sequence of geodesic segments with limit a geodesic line 
7. Now we look at / restricted to 7, this is / o 7 : K — > E such that / o j(t) — > sup / as t — > —00 
and / o j(t) — > —K > —00 as t — > +00. But this is impossible since / o 7 must be strictly concave. 
This proves that / is not bounded below and that / is proper when restricted to T 2 x [0, +00). 

Now we consider C\ = miupxfo} / and C2 < C\ < such that the level set {/ = C2} has at 
least one connected component S in T 2 x (0, +00). Then S is closed and bounded since no sequence 
of points with bounded / can escape to infinity. Therefore S is compact. More explicitly, all level 
sets {/ = C3} with C3 < C2 contained in T 2 x (0, +00) are compact. 

Now we push the level set S to all other level sets by following the flowline <p(x, t) of the vector 
field I grad/l 2 • firstly, the diffeomorphism ip(-, t) brings the level set S C {/ = C2} to the level set 
{/ = C2+ t}, 

f(<p(x, t)) = f(<p(x, 0)) + J ^f(<p(x, «)) ds = f{x) + J (grad /, ±<p(x, a)) ds 

=/w+ r <grad/ 'if^ } ds =^ +t - 

Secondly, the diameter distorsion between these two level sets is bounded. If 7 : [0, 1] — > {0} x T 2 
is a curve on a torus, 

,•,2 r* d 



9 Vt (x)(j,i) = \Mx+ I fo9 v 4x)(%j) ds 





|2 



= \l\x + / ^ g,ad/ (7,7) ds 
JO I grad f\ 2 

Jo I grad/| 2 

Since Hess/ < 0, this implies that ItIL^) < \ j\x so all level sets {/ = C4} with C4 > C2 have 
bounded diameter, and hence are compact and diffeomorphic to S. □ 

Now, the level sets of / are all of them compact and diffeomorphic, thus M = K x {/ = 
c} = K x T 2 and therefore the level sets of / are tori. This allows us to set up a coordinate 
system (r, x, y) £ K x S 1 x S 1 such that the potential function / depends only on the r-coordinate. 
Furthermore, the gradient of / is orthogonal to its level sets, so the metric can be chosen not 
to contain terms on dr ® dx nor dr ® dy. Thus the metric can be written g = u 2 dr 2 + g where 
u = u{r,x, y) and g is a family of metrics on the torus with coordinates (cc, y) parameterized by 
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r. Using isothermal coordinates, every metric on T 2 is (globally) conformally equivalent to the 
euclidean one, thus g = e 2h (dx 2 + dy 2 ) where h = h(r, x, y). These conditions allow us to perform 
computations that reduce to the particular case we studied in Section [2j 

Lemma 9. Consider the metric over IxT 2 

g = u 2 dr 2 + e 2h (dx 2 + dy 2 ) 

where u = u(r,x,y), h = h(r,x,y), and a function f — f(r). Assume that g, f satisfy the soliton 
equation and that g has bounded nonconstant curvature. Then g and f are the ones described 
on the cusped soliton example of Section^ 

Proof. The same riemannian computations as before lead us to the soliton equation 

= Ric + Hess/ + ^g 
= i£ lx dr 2 + ^e 2h E 22 dx 2 + ^e 2h E 33 dy 2 + i£ 12 dr dx + ±E 13 dr dy + \E 23 dx dy 
where 



E n =-u 2 
E 22 = - u 3 

E 33 = - M 3 



fd 2 u 
\dx 2 

( d 2 h 

\dx 2 



dy 2 
d 2 h 

dy 2 



,-2h 



f"u - 2 



Or 



fu - 2 



\ dr J 



(0 + 0)e- 2 ' l -- 2 ( 



dr 



fu 



2( dh 



Or / 



I) 2 a 
dx 2 

d 2 h r} 
dr 2 ( 

d 2 u 

v dy 2 

d 2 h r . 
dr 2 1 



du dh _ 
dy dy 

j du dh 

dr dr ' 

du dh _ 

dx dx 



dh du 
dx dx 



dh du 
dy dy 



or* or 



-2h 



(2f-/0 



-2h, 



du dh 
dr dr ' 



12 


du 

dx 


(dh 
\ dr 


-/')- 


^ dxdr 


13 


du 
~dy 


(dh 
\ dr 


-/')- 


dydr 


23 


du 

~ dy 


dh , 
dx 


du dh 
dx dy 


d 2 u 
dxdy 



E : 

Since the function u never vanishes, nor the exponential does, the soliton equation is the PDE 
system {En — E 22 = E 33 — E\ 2 — E\ 3 — E 23 = 0}. It is convenient to substitute the equations 
E 22 = and E 33 = with the linearly equivalent — \{E 22 + E 33 ) 
E 22 — E 33 = (equation (10)). Then, we get the system 

2 



(5) 

(6) 

(7) 

(8) 

(9) 
(10) 



u 2 e- 2h Au+ -u 3 + f"u~2 



- u 3 e 



-2!i 



du 
dx 
c)u 
dy 



dh 
dr 
Oh 
dr 



r 
r 



du dh 

dy dx 

du dh 

2 2 

dx dx dy dy 



2 

Ah--u 2 e- 2h 
2 

d 2 h 
dxdr 
d 2 h 
dydr 
du dh d 2 u 
dx dy dxdy 

du dh d 2 u 

- 9 1 

dx 2 



dh 
dr 



A' 



dh 
dr 



2 — u 
dr 2 



fu - 2 



du 
dr 

dhY 



(equation (|6j) below) and 
dh 



dr 



r 



d 2 h 



= o, 

du dh 
dr dr 



= 0, 



= 0, 



dy- 



= 0, 
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where A = ^ + ^ is the euclidean laplacian on the zy-surface. We will recover our cusped 
soliton proving that u = 1 and that h(r,x,y) actually does not depend on (x,y). 

We consider first the equations ^ and (10). Since no derivatives on r are present, we can 



consider the problem for r fixed, so u — u(r, •,•) is a function on the xy-torus with metric 
e 2h ( r ''^ [dx 2 + dy 2 ). The function u must have extrema over the torus, since it is compact, so 
there are some critical points (xi,yi) such that |^ |, y , = §f L x y •. = 0. From the equations 

evaluated on a critical point, 3 a * , . = and , . = , , = A ? - so the Hessian 

1 ' axay \(xi,yi) ox- 1 \(xi,yi) "V >{%i,yi) 

matrix (on the rcy-plane) is 



Xi 
Xi 

Suppose that every critical point is nondegenerate, that is, the Hessian matrix is nonsingular with 
Xi 7^ 0. Then the set of critical points is discrete and u is a Morse function for the torus. But 
then the Morse index on every critical point (the number of negative eigenvalues of the Hessian) is 
either or 2, meaning that every critical point is either a minimum or a maximum, never a saddle 
point. Then Morse theory implies that the topology of the xy-surface cannot be a torus (being 
actually a sphere, see [32]). This contradicts that every critical point is nondegenerate, so there is 
some point (xo,yo) such that first and second derivatives vanish. 

We now proceed to derivate the two equations. Equations ^ and ( 10 1 can be written 



(11) Uxy Uyh x ~\- U x hy, 

(12) Hxx ^yy ^lAx^x 2iiyhy , 
using subscripts for denoting partial derivation. Their derivatives are 

"Uxxy — V>xyhx ~t~ V>yh X x ~t~ U X x^y ~t~ V'xhxy t 
V'xyy — Uyyh x ~f~ Uyh X y -\- Uxyhy ~\~ Hxhyyi 
^xxx ^xyy — ^^xx^x ^^x^xx ^U X yhy -\- QtUyh X y , 

Hxxy ^yyy ^^xy^x ^^x^xy ^^yy^y ~t~ ^^y^yy ? 

using the same notation. Evaluated at the point (xq 7 yo), where all first and second order derivatives 
of u vanish, the right-hand side of these equations vanish and therefore all third derivatives vanish. 
Inductivelly, if all n-th order derivatives vanish at (xo,yo)> then the (n — l)-th derivative of the 
equation ( |11[ ) implies that all mixed (n + l)-th order derivatives (derivating at least once in each 



variable) vanish, then the (n — l)-th derivative of the equation ( 12 1 implies that all pure (n + l)-th 
order derivatives (derivating only in one variable) also do; so all derivatives of all orders of u vanish 
at (zojJ/o)- Because u{r,x.y) is a component of a solution of the Ricci flow, it is an analytical 
function (see [5], Ch. 13), so it must be identically constant in (x,y). 

At this point, we can reduce our metric to be g = u(r) 2 dr 2 + e 2h (dx 2 + dy 2 ) with h — h(r, x, y). 
It is just a matter of reparameterizing the variable r to get a new variable, f = J u{r) dr, such that 
u(r) 2 dr 2 = df 2 , so we rename f as r and we can assume that the metric is g = dr 2 + e 2h (dx 2 +dy 2 ) 
with h — h(r, x, y). 

We now look at the equations Q and ^ when u = 1, they imply 

d 2 h d 2 h (( 
dxdr dydr 
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meaning that ^ does not depend on x, y. Finally, looking at equations (|5j), (|6j) when u = 1, we 

Since the left-hand side does not depend on (x,y), nor does the term e~ 2h A h. Recall that a 
two-dimensional metric written as e 2k ( x ' v ^ (dx 2 + dy 2 ) has gaussian curvature K = —e~ 2h A h. So 
the xy-tori have each one constant curvature, and the only admitted one for a torus is K = 0. 
Hence h only depends on r and the equations turn into the system ^ , that we already studied for 
the example of the cusp soliton. The rest of the uniqueness follows from the discussion on Section 

m □ 

4. Evolution of curvature 

On this last section we expose the property anounced in Theorem |1.3[ derivated from the 
opposite effects of the diffeomorphism and the homothety for the evolution of the metric. Recall 
that (M,g(t)) is the (soliton) Ricci flow defined on M — R x T 2 and for t € (—1, +oo), such that 
<?(0) = go where go is the metric constructed in Theorem |l.l[ and let us denote R = R(t) the scalar 
curvature of g(t). We want to show that the growth of the curvature along (M,g(t)) changes sign 
for values of t far enough of —1, but is positive everywhere along the manifold for values of t close 
enough to —1. 



Proof (of Theorem 1.3). The evolution of the soliton metric under the Ricci flow is 

ff(t) = (* + iKW 

where go is the metric constructed in Theorem 1 1 . 1 1 and <j>t(ro, Xq, yo) = (r(t),x ,y ) with 

f r(t) = f'(r(t)) = F(r(t)) 
1 r(0) = r . 

Then, since R[g ] = —AH' - QH 2 (by Lemma[l|, 

1 „ r „ 1 



2\ 



SO 



jRW)](r oM = Jj^y>i-4H' 6H 2 ) + j^(-AH"F 12HH' F)\ r=r(t) 

1 [{2HF - H 2 + 1) + (t + 1)F 2 (-2HF + 2H 2 - 1)] | = 



(t + l) 2 LV ' V ' V ' l \v=r(ty 

Thus, the zeroset {^R[g{t)\ = 0} defines, for each t, an algebraic curve on the ffF-plane. If 
the solution curve S intersects this zeroset curve, then the soliton changes the growth sign of the 
curvature at some point. Otherwise R is everywhere monotone. 

For the rest of the proof, we rename the variables so our system Q is 



(13) 



x = xy — 2x 2 + | 
y = 2xy - 2x 2 + \ 



where x = x(r), y = y(r), the curve S is the separatrix solution of the system ( 13 1 emanating from 
the critical point (^,0) towards the vertical asymptote x = 0, and 

(14) C t = (2xy - x 2 + 1) + (t + l)y 2 (-2xy + 2x 2 - 1). 
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The question is whether {Ct = 0} intersects S, for each t £ (—1, +00). Figure [5(a) represents the 



solution curve S together with the curve Ct for t = 10. This gives some evidence that for big values 
of t there is an intersection point of the curves, but it is not clear for small or negative values of 
t. The issue is that C t has an asymptote and whether it approaches the infinity at the right hand 
side or at the left hand side of the curve S. In order to study these guesses, we perform again a 
projective change of variables, equivalent to assume that our phase portrait lies on the {z = 1} 
plane (with coordinates (x, y, 1)) of the xyz-space, and we project perspectively from the origin to 
the {y = —1} plane (with coordinates (x, —l,y)). 

x = — - J x = I 

v = -I ' \ y = -I • 

This change of coordinates has the effect of bringing the point at the infinity on the vertical 
asymptote to the new origin of coordinates, the old line at infinity to the horizontal axis, and the 
old (projective) line y = to the new line of the infinity. We won't keep track of the tilde notation 
and use again x, y as coordinates. 

After this change, the system turns into 



f x. _ -8a -ix s +xy -2x+y 

(15) \y = -2x- 2x 2 \ \y 2 

that is equivalent (has the same orbits) to the system 

(16) 

and the curve Ct turns into 



x = —4a;" — ix" + %xy — x + \y 2 

2„, 1 1„,3 



y = -2xy - 2x l y + ±y 4 
-2xy 2 - x 2 y 2 + y 4 + (t + l)(2x + 2x 2 - y 2 ) 



y 4 

that has the same zeroset as 

(17) C t = -2xy 2 - x 2 y 2 + y 4 + (t + l)(2x + 2x 2 - y 2 ). 

See Figure |5(b)| In particular, we can check that Ct now passes through the origin, and this 
confirms that the original Ct in (14) had a vertical asymptote. Nevertheless, it is not yet clear 
which curve lies at which side near the contact point. In order to investigate this behaviour, we 
perform some algebraic blow-ups at the contact point. Recall that an algebraic blow-up is a change 
of variables from the old (a;, y) to the new (x , y) given by 



x = xy 



y = y l v = y 

The mapping tp : (x,y) 1— ^ (x,y) is a birrational map, which restricts to a diffeomorphism in all 
points except at (x,y) = (0,0), where </3 -1 ((0,0)) = {y = 0} is a (projective) line called the 
exceptional divisor of the blow-up. The exceptional divisor is in correspondence with the space 
of directions of the old origin, thus we "pick out a point and substitute it with a projective line" . 
Again, we won't keep track of the tildes. Two curves intersecting with normal crossing at the 
origin are transformed in this way to two separated curves; two curves tangent at the origin, 
when transformed, still intersect, but their contact order is decreased. Since all the curves we are 
involved with are analytical, after a finite number of blow-ups the process finishes separating the 
curves. With the exception of the multiply blown- up line {y = 0}, all the remaining phase portrait 
is diffcomorphic to the original one, so any intersecting point other than the origin will still be 
present in the blown- up portrait. 
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The process can be algorithmically carried on. We consider the vector field of the system ( 16 1 



The solution S intersects the {y = 0} axis at the (one) critical point of the vector field, that can 



be symbolically computed. We consider also the curve Ct in (17) after the chart change. This 
is a polynomial in x, y and its intersection with {y = 0} can also be computed symbolically. We 
perform the change of variables corresponding to the blow-up, and ocasionally translate the new 
intersection point to the origin again. We compute both intersection points with {y = 0} and 
iterate up to when the two results disagree. Let us remark that this process can be carried out 
by a symbolic algorithm, so there is no numerical approximation involved. See Figure [5] for a 
numerical visualization. 

Once this is done, we find that after six blow-ups the critical point of the system is located at 
(0,0), and the intersection of C t with the {y = 0} line is at (|^rf,0). Thus, generically the curve 
Ct and the solution curve S have order of contact five at the infinity in the original phase portrait 



of (13). In the case t = 0, both points of intersection agree, so further blow-ups are needed. It 
turns out that the tenth blow-up separates the points, and when locating the critical point at (0, 0), 
the intersection of C t with {y = 0} is at (|, 0). Thus the curve Cq has order of contact nine with 



the solution S at the infinity in the original phase portrait of (13). 

Now we recall a couple of properties of the blow-ups: firstly, curves crossing the origin with a 
slope A are blown-up to curves crossing the {y = 0} line at x = -r, in particular positive slopes 
are sent to the x > half-line, and negative slopes to the x < half-line. Secondly, the blow-up 
preserves orientation of horizontal lines in the upper half-plane, and reverses it in the lower one. 

The change of charts we performed before the blow-ups preserves the orientation of all horizontal 
lines, but exchanges the lower and the upper half-planes. In summary, the relative position (left 



and right) of Ct and S on the lower half-plane of the original phase portrait of (13), is the same 
as on the upper half-plane on the phase portrait after all the blow-ups. 

Therefore we can deduce that for t £ (—1,0) the curve C t approaches the infinity at the asymp- 
tote from the left-hand side of the separatrix S, and for t £ [0, +oo) it approaches from the right. 
Given that this component of the Ct curve has always points at the left-hand side of S, we can 
deduce that Ct intersects S at least at one point (other than the infinity) for t > 0. 

Furthermore, this finite intersection point depends continuously on t, so a small perturbation 
on t will still make the two curves intersect. Thus, there exist a small 6 < such that for t > 5 
the function C t changes sign along the separatrix S. Actually, for 8 < t < the sign must change 
at least twice. 

Finally, we show that this change of sign of C t does not happen for some t close enough to 
— 1. This is due to the fact that for such t the curve {Ct = 0} is a barrier for the separatrix. We 



compute the normal vector to the curve Ct in (14), 



grad C t = (2y-2x+(t + l)y 2 (-2y + 4:x),2x + 2(t + l)y(-3xy + 2x 2 - 1)) 



and compare with the vector field of the system ( 13 ) 



V = [xy - 2x 2 + X -,2xy - 2x 2 + X - 

Their scalar product is 
(grad C u V) = 

= -y (-2xy + 2x 2 - 1 + (t + l)(2xy 3 + 4x 2 y 2 + y 2 ~ 12x 3 y + bxy + 8x 4 - 6x 2 + 1)) 



Restricted to the curve {Ct = 0}, this simplifies by substracting the equation yCt = 



(18) M X ,V) : = (grad C t , | c = -y (x 2 + (t + l)(6x 2 y 2 - 12x 3 y + 5xy + 8x 4 - 6x 2 + 1)) 
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Fortunately Ct — is a second order equation for x, thus it is easy to select the apropriate isolation 
x — x{y) corresponding to the branch on x > 0, y < 0, 



x = x{y) 



-y + (t + 1) y 3 + ^jy 2 -2 (t + l)y 4 + {t + if if - 3 (< + 1) y 2 + 2 (t + I) 2 j/ 4 + 1 

2 (i + l)y 2 -l : 



and substitute it on (18) (although unfortunately, the explicit expression is quite ugly), 

*t(y) := ipt{x(y),y) = ... 



This is just a real function ^ t 



—oo, 



(the upper bound on the domain is the 



negative solution of Ct(0,y) = 0). This function gives for each value of y the scalar product 
between the normal vector to the C t curve (pointing rightwards) and the vector field of the system 
at the point (x(y),y). If this function is strictly positive for values of t close to —1, this implies 
that C t is a barrier for the separatrix. For, the separatrix S emanates from (1/2,0), which is on 
the right hand side of C t , and if S touched C t , then its tangent vector (the vector field of the 
system) would be pointing to the same region separated by C t . Otherwise, if the function ty t fails 
to be positive, then Ct fails to be a barrier. As shown in Figure |4j Ct is a barrier for t = —0.7 but 
it is not for t = —0.2. The given value of —0.7 is just an example, it could be checked for smaller 




(a) t - 



-0.7 



(b) t = 



-0.2 



Figure 4. Graph of ^t{y) for two values of t. 

values. However, it is not immediate to tell which is the critical value, since the failure of Ct to be 
a barrier does not ensure an actual crossing of Ct and the separatrix S. 

We can therefore say that the scalar curvature is always negative, and that (for instance) for 
t = —0.7 the soliton is evolving with the scalar curvature everywhere increasing, whereas for t > 
the soliton has regions where the scalar curvature is increasing and regions where it is decreasing. 
Any fixed point eventually belongs to the region of increasing curvature and therefore any point 
eventually tends to zero curvature, due to the dominance of the expanding effect. □ 
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Figure 5. Analysis of the tangency at infinity of the separatrix S (bold line) and 
the curve {C t — 0} with t = 10. Shaded areas correspond to C t > 0. 
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